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Abstract 



^ ■ We analyze the supersymmetric (SUSY) QCD contribution to the h'^bb 

coupling at one loop in the Minimal Supersymmetric Model (MSSM) in the 
decoupling limit. Analytic expressions in the large SUSY mass region are 
r> I derived and the decoupling behavior of the corrections is examined in various 

d ' limiting cases, where some or all of the SUSY mass parameters become large. 

We show that in the decoupling limit of large SUSY mass parameters and 
large CP-odd Higgs mass, the h'^bb coupling approaches its Standard Model 
value at one loop. However, the onset of decoupling is delayed when tan/? 
is large. In addition, the one-loop SUSY-QCD corrections decouple if the 
masses of either the bottom squarks or the gluino are separately taken large; 
although the approach to decoupling is significantly slower in the latter case. 
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1 Introduction 

Once a light CP-even Higgs boson is discovered, it will be crucial to measure 
as many of its couplings as possible with the highest accuracy possible. By 
measuring the Higgs couplings to gauge bosons, one can learn whether only 
one Higgs boson is involved in electroweak symmetry breaking dynamics. 
Moreover, the Higgs couplings to vector bosons are sensitive to the possible 
existence of non-doublet isospin structure in the Higgs sector. By measuring 
the Higgs couplings to fermion pairs, one can learn whether the Higgs boson 
is responsible for fermion mass generation. Knowledge of the trilinear and 
quartic Higgs self-couplings, although extremely difficult to obtain, would 
allow one to reconstruct the Higgs potential and directly test the mechanism 
of electroweak symmetry breaking. Finally, if the couplings can be measured 
at the level of the radiative corrections, one could then derive significant 
constraints on new physics beyond the reach of the present accelerators. A 
detailed study of radiative corrections to the Higgs couplings would be es- 
pecially important if a light Higgs boson were discovered in the mass range 
predicted by the minimal supersymmetric extension of the Standard Model 
(MSSM), but supersymmetric (SUSY) particles were not found. In this case, 
the precise experimental determination of Higgs couplings could provide in- 
direct information about the preferred region of SUSY parameter space. For 
example, one could predict (in the context of the MSSM) whether the data 
favored a SUSY spectrum below the 1 TeV energy scale. 

It is well known that the tree-level couplings of the lightest MSSM Higgs 
boson (/i^) to fermion pairs and gauge bosons tend to their Standard Model 
(SM) values in the decoupling limit. Ma ^ Mz Q, where Ma is the mass of 
the CP-odd neutral Higgs boson (A^) of the MSSM. As a consequence of this 
decoupling, distinguishing the lightest MSSM Higgs boson in the large Ma 
limit from the Higgs boson of the Standard Model (SM) will be very diffi- 
cult. Formally, the decoupling of all SUSY particles (including the radiative 
corrections) implies that in the effective low-energy theory, all observables 
tend to their SM values in the limit of large SUSY masses and large Ma- 
It has been shown that all of the SUSY particles in the MSSM, including 
the heavy Higgs bosons H^, ^4° and i7*, decouple at one- loop order from 
the low-energy electroweak gauge boson physics p. In particular, the con- 
tributions of the SUSY particles to low-energy processes either fall as inverse 
powers of the SUSY mass parameters or can be absorbed into counterterms 
for the tree- level couplings of the low-energy theory 0. As a result, the 
radiative corrections involving SUSY particles go to zero in the asymptotic 
large SUSY mass limit. Our aim is to determine the nature of the decoupling 
limit at one-loop for the couplings of h^ to SM particles. 

In this paper, we focus on the h^ coupling to bb. This coupling determines 
the partial width of h^ -^ bb, which is by far the dominant decay mode of 
h^ in most of the MSSM parameter space. Because this decay is dominant, 
accurate knowledge of the h^bb coupling is very important for Higgs boson 



searches. At LEP and the Tevatron, the primary Higgs search channel is 
h^ -^ bb. The experimental reach of Higgs boson searches at the upcoming 
Tevatron Run 2 depends critically on the h^ -^ bb branching ratio |^. In 
contrast, the Higgs boson searches at LEP do not depend as critically on the 
/i° — >• bb decay. At LEP, there is sufficient cross-section to detect the Higgs 
boson in multiple channels. Moreover, even without observing the Higgs 
decay products, the Higgs boson mass can be reconstructed by detecting the 
recoiling Z boson in e~^e~ — > Zh^. At the CERN LHC, the primary Higgs 
search channel in the mass region below 130 GeV is the rare decay /i° — > 77. 
The Higgs event rate in this channel is affected strongly if the total width 
of h^ is modified due to corrections to the dominant bb decay mode [^. The 
same holds true for other search channels at the LHC such as h^ -^ r^r^ 0. 

In this paper we study the MSSM radiative corrections to the h^bb cou- 
pling at one loop, to leading order in a^, and we analyze in detail their 
behavior in the decoupling limit. [| These corrections are due to the SUSY- 
QCD (SQCD) sector and arise from gluino and bottom-squark (sbottom) 
exchange. Because of the dependence on the strong coupling constant, these 
are expected to be the most significant one-loop MSSM contributions over 
much of the MSSM parameter space. Potentially significant contributions can 
also arise from the SUSY-electroweak sector (the most significant of which 
are proportional to the Higgs-top quark Yukawa coupling); we will address 
these corrections elsewhere and do not consider them here. 

The SQCD corrections to the h^bb coupling were ffist calculated using 
a diagrammatic approach in ref. (which also contains results for the 
SUSY-electroweak corrections). These corrections have also been obtained 
in refs. 0,[l3. The SQCD corrections were also computed in an effective La- 
grangian approach in ref. [Q, using the SUSY contributions to the 6-quark 



self energy obtained in refs. |T2,|T^ and neglecting terms suppressed by in- 
verse powers of SUSY masses. 

The radiatively-corrected h^bb coupling depends on the CP-even Higgs 
mixing angle a. At tree- level, this mixing angle is determined by fixing tan/3 
and Ma- At one-loop order, there are no 0{as) corrections to this mixing 
angle. As a result, working to leading order in a^, we may employ tree- level 
relations for a in our computation of the h^bb coupling. This procedure is 
no longer adequate once one-loop SUSY-electroweak effects are included. In 
the latter case, the one-loop radiative corrections to a must be taken into 
account, as described in refs. [§|,|^,|13|. These papers show that the interplay 
between the radiative corrections to the mixing angle and to the h^bb coupling 
can be very important for Higgs collider phenomenology, particularly in the 
case when the branching ratio for h^ -^ bb is suppressed.^ This is most 
easily seen as follows. When radiative corrections to the mixing angle a 



^The radiatively-corrected h'^bb coupling in the decoupling limit has also been consid- 
ered previously in refs. |Q,0,Ql- 

^The relevance of the suppressed h'^bb coupling for phenomenology has also been em- 
phasized in refs. |,|l|,|l|]. 



are included, it becomes possible to tune this angle to zero independent 
of the value of tan/? by varying the SUSY parameters. At a = 0, the 
tree-level couplings of h^ to bb and t~^t~ vanish, as do the ordinary QCD 
corrections |jl6| to the h^bb coupling. However, because the SQCD corrections 



to the h^bb coupling include contributions from diagrams involving the /i° 
coupling to sbottoms, they remain nonzero at a = 0. As a result, the h^T^T~ 
coupling goes to zero at a different point in SUSY parameter space than the 
h'^bb coupling does [|, |],|13|. We will come back to these issues and study 
the approach to decoupling of the SUSY-electroweak corrections in a later 
paper. 

In some regions of the MSSM parameter space, the SQCD corrections to 
the h^bb coupling become so large that it is important to take into account 



higher-order corrections. This has been carried out in refs. |17, 18| by re- 
summing the leading tan j3 contributions to all orders of perturbation theory 
using an effective Lagrangian approach. This resummation is not important 
in our present work because we are interested in the decoupling limit in which 
the one-loop corrections to the h^bb coupling are small. 

In this paper we analyze the full diagrammatic formulae for the on-shell 
one-loop SQCD corrections to the h^bb coupling. We perform expansions in 
inverse powers of SUSY masses in order to examine the decoupling behavior 
when the SUSY masses are large compared to Mz- The SQCD corrections 
depend on a number of different SUSY mass parameters, and the relative 
sizes of these masses affect the manifestation of the decoupling. To remain 
as model-independent as possible, we make no assumptions about relations 
among the SUSY parameters that may arise from grand unification or specific 
SUSY-breaking scenarios. We consider the soft-SUSY-breaking parameters 
and the /i parameter to be independent parameters whose magnitudes are 
all of order 1 TeV. 

In this paper, we demonstrate that in the limit of large Ma (in this 
limit one also has MHO,Mfj± ^ Mz) and large sbottom and gluino masses 
{Mf^,,Mg ^ Mz), the SM expression for the h^bb one-loop coupling is re- 
covered. That is, the SQCD corrections to the h^bb coupling decouple in 
the limit of large SUSY masses and large M^. In particular, we examine 
the case of large tan/?, for which the SQCD corrections are enhanced. This 
enhancement can delay the onset of decoupling and give rise to a significant 
one-loop correction, even for moderate to large values of the SUSY masses. 

This paper is organized as follows. In Section ^ we define our notation 
and briefly review the Higgs and sbottom sectors of the MSSM. In Section | 
we give the exact one loop formula for the SQCD corrections to the h^bb cou- 
pling. In Section ^ we derive analytic expressions for the SQCD corrections 
in the hmit of large SUSY masses. We analyze the decoupling of the SQCD 
corrections for various hierarchies of mass parameters, and numerically com- 
pare the analytic approximations to the exact one- loop result. In Section ^ 
we summarize our conclusions. Finally, the Appendix contains expansions of 
the one-loop integrals used in our calculations. 



2 Higgs and sbottom masses in the MSSM 

In the MSSM, the parameters of the Higgs sector are constrained at tree-level 
in such a way that the Higgs masses and mixing angles depend on only two 
unknown parameters. These are commonly chosen to be the mass of the CP- 
odd neutral Higgs boson A^ and the ratio of the vacuum expectation values 
(vevs) of the two Higgs doublets, tan/5 = f2/fi- (For a review of the MSSM 
Higgs sector, see |^.) In terms of these parameters, the mass of the charged 
Higgs boson H"^ at tree level is Mjj± = M\ + M^ , and the masses of the 
CP-even neutral Higgs bosons h^ and H^ are obtained by diagonalizing the 
tree-level mass-squared matrix, 

2^ / Ml sin2/3 + M|cos2/? -(Ml + M|) sin /3 cos /3\ 

1^ - (Ml + M|) sin /? cos /3 Ml cos^ /3 + M| sin^ /3 ) ' ^^'^^ 



The eigenvalues of this matrix are, 
1 



M 



HO,hO 



2^2 



M% + M'z±^ {Ml + MlY - AMlMl cos2 2/5 



(2.2) 



with M/jO < Mfjo. At tree- level, M^o < Mz\ cos 2/3 1; this bound is saturated 
at large Ma- We choose a convention where the vevs are positive so that 
< /5 < 7r/2. The mixing angle that diagonalizes A^^ is given at tree- level 
by 



M\ + Ml 
Ml -Ml 



tan 2a = tan 2/3 —I -| . (2.3) 



In the conventions employed here, — 7r/2 < a < (see ref. [^ for further 
details). From the above results it is easy to obtain: 

oosHP - a) = ^^M^^tl^^ (2 4) 

'°'^^ "^ MKM^o-M^o)- ^'-^^ 

In the limit of Ma ^ Mz, the expressions for the Higgs masses and mixing 
angle simplify and one finds 

Mlo ~ M^ cos^ 2/3 , 
Mjjo ~ Ml + M| sin^ 2/3 , 

cos'if3-a)^ ^^^, ^ . (2.5) 

Two consequences are immediately apparent. First, M^ ~ Mj^o ~ M/^±, up 
to corrections of O^M^/Ma)- Second, cos(/3 — a) = up to corrections of 
0{Mz/Ma)- This limit is known as the decoupling limit because when Ma 
is large, one can define an effective low-energy theory below the scale of Ma 
in which the effective Higgs sector consists only of one light CP-even Higgs 



boson, h^, whose couplings to Standard Model particles are indistinguishable 
from those of the SM Higgs boson [|l| . From eq. p.5|, one can easily derive: 



2M^ ( M'^\ 

cota = -tan/?-— ^tan/?cos2/5 + C — I . (2.6) 

When radiative corrections to the CP-even Higgs mass-squared matrix 
are taken into account, the upper bound on My,o increases substantially to 
Mf^o < 135 GeV (assuming all supersymmetric particles are no heavier than 
about 1 TeV), and corrections to a become substantial for low Ma- These 



corrections are well known |12,21-26| and the leading contributions have 



been computed up to two-loop order. In this paper we consider only the 
contributions to the hPbb coupling of order a^ at one loop. Because the 0{as) 
contributions to the CP-even Higgs mass-squared matrix only first arise at 
the two-loop level, the radiative corrections to this matrix are irrelevant to 
our present work. (In contrast, they do contribute to the one-loop SUSY- 
electroweak corrections to the h%h coupling.) 

From direct searches at LEP the MSSM h^ and A^ masses are constrained 



to be A4o > 88.3 GeV and Ma > 88.4 GeV ^. For a range of values of 
tan/5 close to one, the theoretical upper bound on M^o is lower than the 
experimental lower bound, so the corresponding region of tan (3 can be ruled 
out. Because of the radiative corrections, the variation of the upper bound 
depends primarily on the precise value of the top quark mass and the mixing 
in the stop sector. For the conservative choice oirrit < 179 A GeV and mixing 
in the stop sector that maximizes the upper bound on M^o, values of tan/5 
between 0.8 and 1.5 are excluded [ p7| . 

We now discuss the parameters of the sbottom sector. The tree-level 
sbottom squared-mass matrix is: 

^f - ( »1 "t ) ■ <"' 

where we use the notation, 

Xh = Aj, — fitanP , 

Ml = Ml + ml + Ml{ll - Q.sl^) cos2/3, 
Ml = Ml + ml + MlQbslr cos 2/3 . 



(2.8) 



Here /| = —1/2 and Qh = —1/3 are the isospin and electric charge of the 
6-quark, respectively and sw = sin6'vy. The parameters Mq and M^ are the 
soft-SUSY-breaking masses for the third-generation SU(2) squark doublet 
Q = (Jl, bi) and the singlet D = 5|j, respectively. Af, is a soft-SUSY-breaking 
trilinear coupling and /x is the bihnear coupling of the two Higgs doublet 
superfields. The sbottom mass eigenstates are 

bi = cos 6ibL + sin O^bji; b2 = -sin O^bL + cos Oibn, (2.9) 



where — 7r/4 < 9^ < 7r/4 is defined so that 61 (62) is predominantly bi (&r) 
The sbottom mass eigenvalues are then given by 



Mi 



where0 



Ml + M'^± cTLRJiMl - MlY + AmlXl 



^LR = ^m{Ml - Ml) 



(2.10) 
(2.11) 



and the mixing angle 9i is given by 



cos 29r 



sin 297 



Ml -Ml 
Ml -Ml 

2mbXb 
M? -Ml 



Note that in our conventions, 



Ml, > 



the sbottom masses switches if cr^^ < 0. 



Mi, if a^H 



(2.12) 



> 0, whereas the order of 
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the sbottoms must be heavier 
than about 140 GeV, assuming that the mass of the lightest neutralino x5 
is less than half the mass of the lighter sbottom. For heavier neutralinos, 
the Tevatron searches lose efficiency. In this region the direct searches at 
LEP 1^ place a lower bound on the sbottom masses of about 85 GeV. The 
LEP bounds are valid only for b — Xi niass splittings larger than about 5 
GeV, so that the decay mode b — > 6x1 is kinematically accessible. 

The limits on the gluino mass Mg are more model-dependent. If one 
assumes relations between the gaugino masses such that they unify at the 
GUT scale, then Mg is constrained from direct searches at the Tevatron to 



be greater than 173 GeV, independent of the squark masses ^0 



3 SQCD corrections to h^ -^ bb 

The h'^bb coupling is given at one-loop level to order a^ by 

ghbb = ghbb + ^Qhbb + ^9hbb — 9hbb (1 + Aqcd + ^sqcd) 



(3.1) 



where gubb is the one-loop coupling, ghbb is the tree-level coupling, ^g^b^ i^ 



the radiative correction from pure QCD |[TO], and ^g^^b is the one-loop 
SQCD contribution. 

The tree-level h^hb coupling is given by 



Qhbb 



grrib sm a 
2Mwcos(3 



(3.2) 



■'if Ml = Mr, then cr^^ is not well-defined. In the present context, a useful convention 
is to set aj^j^ = a^ [where a^ = sgn(X;,)] if Ml — Mr. Nevertheless, one can check 
that our final expressions for the radiative corrections in Section are independent of this 
choice. 



h^ 






hP 




h' 




h'- — 




Figure 1: Feynman diagrams for the SQCD corrections to the h^bb couphnj 
The vertex marked with the X refers to the one-loop h^bb counterterm. 



Note that in the hmit of large Ma, sin a — > — cos/3 and ghbb tends to the 
SM coupling, (yff^^ = —gmb/{2Mw)- The one-loop corrections to the h^bb 
coupling modify the h^ — > bb decay width as follows, keeping only correction 
terms of 0{as): 



T{h^ -^ bb) = r(/i° -^ 66) (1 + 2Aqcd + 2A 



SQCD) 



(3.3) 



where T is the one-loop partial width and F is the tree-level partial width. 

The SQCD contribution to the h^bb coupling comes from diagrams in- 
volving the exchange of virtual gluinos (g) and sbottoms (bi), as shown in 
fig. |I|. We have 



^gSQCO 



ihbb 



{S ghbb) 3^^^ 



+ (S ghbb) 2'^'^'^ 



+ iS9hbb)x^^^ 



(3.4) 



consisting of contributions from the vertex correction, the 6-quark wave func- 
tion renormalization, and the counterterm from the renormalization of the 
6-quark Yukawa couphng, respectively. To compute the one-loop Yukawa 
counterterm contribution, we note that the Higgs wave function, the vevs 
(and hence tan/5) and the parameters g, Mw and a receive no 0{as) cor- 
rections at one- loop. Thus, to leading order in a^, {Sghbb)x can be easily 



obtained from eq. p.2| and depends only on the 6-quark mass counterterm as 
follows: 



{^9hbb)^^^ 



Qhbb 



{6mt 



\SQCD 



nib 



(3.5) 



In eq. |3.5| , [Smf,)^'^'^^ is the SQCD contribution to the 6-quark mass coun- 
terterm, which is fixed by defining rrih to be the pole of the one-loop 0{as) 
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6-quark propagator. This is the on-shell renormahzation scheme. 

We have computed the various contributions to Sgf^jff^ [see eq. ^^ . The 
contribution of the one-loop vertex is given by: 



{5ghbb)l^^^ 



Qhbb 



37r' 



2M7COs/?sin(Q; + /5) , , a ^ 9 x 



rrib 



sma 



X [rribCn + MgS2bCo] (ml Ml,, ml M?, Ml , M^J 



+ 



2M| cos/3sin(a + /3) ,^ft n ^ ? x „ 

£ t^_A ^{I'^sl + QbslrC2b) + 2mh - YbS2b 

rrib sm a 

X Kdi - MgS2bCo] K, M^o, mlM? M? M-^ 



62 ' b2> 



rrib sm a 



X 



2M5C2,Co(m2, Ml,,ml- Ml Ml, Ml) 



(3.6) 



where Cb = cos 6*^, C2b = cos26'^, Sb = sin 6*^, etc., and Yb arises in the Higgs 
couphng to sbottoms: 

Yb = Ah + fi cot a . (3.7) 

The contribution from the 6-quark self-energy and the h^bb vertex coun- 
terterm is given by 



■.SQCD 



{^9hbb)2 + {^9hbb)x 



SQCD 



Qhbb 

37r I nib 



a 



S2b 



Bo{ml;MlMl)-Bo{ml;MlMl 



-2ml 



bi ^^-^gi ^^%2' 



B[{mi;Ml,Ml) + B[{mi-Ml,Ml) 
2mbMgS2b 



B',{ml-MlMl)-B',{ml;MlMl: 



(3.8) 



Our notation for the loop integrals Bq, B'q, B[, Cq and Cu is defined in 
the Appendix. We have checked that our results are in agreement with the 
calculations of ref . [H . 



4 Analytic and numerical results 

We now analyze the decoupling behavior of the SQCD corrections to the 
h^bb coupling. We derive approximate analytic expressions for the SQCD 
corrections in the limit of large SUSY mass parameters and explore the nature 
of the decoupling limit. 



We define our expansion parameters as follows. Since we are interested 
in the limit of large SUSY mass parameters, we consider all the soft-SUSY- 
breaking mass parameters and the fi parameter to be of the same order (col- 
lectively denoted by Msusy) and much heavier than the electroweak scale. 
That is, 

Msusy -- Ml -^ Mr^ Mg r^ ^ r^ A^:^ Mew, (4.1) 



where Mi and Mr are defined in eq. p.8| . We give expansions of the SQCD 
corrections to the h^hb coupling in inverse powers of the SUSY mass pa- 
rameters, up to order M|;^y/M|j/_5y. We consider Mz-, Mho, rubtanP, and 
rubcota to all be of order Mew- We neglect small contributions of order 
ml/MgjjgY aiid rribMEw /MgjjgY that are not enhanced by tan/? or cot a. 
The expansions of the loop integrals are given in the Appendix. There are 
two possible extreme configurations of the sbottom mass-squared matrix that 
we must separately consider: maximal and near-zero mixing. 

Maximal mixing {Oi ~ ±7r/4) between hi and h^ arises when the splitting 
between the diagonal elements of the mass-squared matrix is small compared 
to the off-diagonal elements: |M| — M^| -C mb\Xb\. Because of the tan/? 
enhancement in X^, m^Xb is of order MewMsusy- In this case we consider 
|M| — M|,| to be of order M'^y^, so that the mass splitting between the two 
sbottoms is small compared to their masses and we must take care to treat 
it properly in the expansions. We consider this case in Section |4.1| . 

Near-zero mixing between hi and hR arises when the splitting between the 
diagonal elements of the mass-squared matrix is large compared to the off- 
diagonal elements, |M| — Mr\ 3> mb|Xfe|. This is the case usually considered 
in the literature, because M^ and Mr depend on two different soft-SUSY- 
breaking parameters Mq and M^, respectively, and the 6-quark mass in the 
off-diagonal elements is small. In this case the mass splitting between the 
two sbottoms will be of the same order as their masses {i.e., {Ml — M|,| is 
of order Mgrjay) and this has to be treated properly in the expansions. We 



consider this case in Section 4.2 



4.1 Maximal bi — hR mixing 

Maximal mixing in the sbottom sector arises when |M| — Mr\ <^ mb\Xb\. In 
this limit, we can expand the sbottom mass-squared eigenvalues in powers of 
the small parameter (M| — Mr) / mbXb (which is of order Mew /Msusy) as 
follows: 

Ml^^ ^ M| ± A^ , (4.2) 

where we have defined 

Ml = \{Ml + Ml) = \{Ml+Ml) 
, (M|-M|)2 



A = a^p,mb\Xb 



SmlXi 
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(4.3) 



Here Mg is of order Mg^gY while A^ is of order MewMsusy- Expanding 
the expressions for the mixing angle in terms of the same small parameter, 
we obtain 



cos 26r 



Ml -Ml 



2mhXi 



b^b 



sin2%~a^^a^ 



1- 



(M 



Mir 



8mlX^ 



(4.4) 



where ax = sgn(Xfe). Expanding eqs. |3]^ and ^]8| to order Ml^r/MgJJgYy we 
find 



A 



a. 



SQCD 



rml tan^ (3 

2MJ 



-I^Mg 

Ml 

cot a 
tan/9 



(tan P + cot a) /i (R) 



YbMgMlo 
12Mi 



MR) 



/4(^)-^fn-2^^V3(^) 



Ml 



+ - 



2 M| cos /5 sin (a + /5) . 



3 Ml 



sma 



II MR) + 



MgXk 

Ml 



tan/5 



f2{R)]+0 



'rUbM, 



EW 



M^ 



,(4.5) 



where R = Mg/Ms- The functions fi{R) arise from the expansions of the 
loop integrals and are given in the Appendix. They are normalized so that 
/j(l) = 1. Note that terms of order (M| — M|,)^/(m^X^) cancel exactly in 
the leading order of the large Msusy expansion [eq. |4.5|| . 

The first term in eq. ff3| is zeroth order in Msusy- That is, if the ratios 



between the SUSY parameters are fixed and the SUSY mass scale is taken 
arbitrarily heavy, this term remains constant. This non-decoupling behavior 
has been pointed out previously in refs. @, Q. If the SUSY mass scale is 
much larger than M^, then one may define a low-energy effective theory by 
integrating out the SUSY particles. This low-energy effective theory contains 
two Higgs doublets, whose couplings to fermions are unrestricted {i.e., each 
Higgs doublet couples to 6ot/i up-type and down- type quarks), characteristic 
of the so-called general type-III model [^ . 

The remaining terms are of order Mly^/Mlugy- I^ contrast to the first 
term, they depend on A^ (through X^ and Yj,). However, the contribution 
proportional to A^ is not enhanced when tan/5 (or cot a) is large, and so is less 
significant at large tan/5 than the contribution proportional to /i. Neglecting 
all terms that are not enhanced by large tan/5 or cot a, we find that ^sqcd 
is proportional to the product iiMg. Because of this, for large tan/5 the sign 
of AsQCD can be used as a test of the anomaly-mediated SUSY breaking 
32| , which predicts a negative Mg |]33|. Of course, the sign of n 



scenario 



must be determined from another SUSY process for the sign of Mg to be 
extracted. 
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4.2 Near-zero bi — Br mixing 



Near- zero mixing in the sbottoni sector arises when |M| — M^| 3> mftlXftl. 
This corresponds to taking the difference between the physical sbottoni masses 
to be of the same order as the masses themselves. In this case we write 
our results in terms of the physical sbottom masses and expand in powers 
of the small parameter mbXb/{M^ — M^ ), which we take to be of order 



62^ 



Mew/Msusy- The mixing angle is then given by eq. |2.12| , from which one 
easily derives the expansion 



cos 2^^ ~ 1 



2mlX^ 
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b2' 



(4.6) 



Expanding eqs. ^]6|and|3]^ to order M^y^/MgjjgYi ^'^^ writing the result 
in terms of the physical sbottom masses, we find: 
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(4.7) 



where Ri = Mg/M^^ {i = 1,2). The functions /ij(i?i,i?2) and /i,5(-Ri) arise 
from the expansions of the loop integrals and are given in the Appendix. 

As in the case of maximal sbottom mixing, the first term in eq. ^^ is 
zeroth order in Msusy- The remaining terms are of order M'^^/M'gjjgY- ^^ 
in the previous section, if we neglect all terms that are not enhanced by large 
tan/3 or cot a, we find that the dependence on A^, drops out and Asqcd is 
again proportional to the product fxMg. 



4.3 The approach to the decoupUng hmit 



If we take all SUSY mass parameters large at fixed Ma in eqs. ^^ and |477| , 
then Asqcd tends to a nonzero constant; i.e., the SQCD corrections do 
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not decouple. However, we are especially interested in the case where both 
MsusY and Ma are large. In eqs. |4.5| and |4.7|, the terms of zeroth order in 



MsusY are proportional to tan/? + cot a. From eq. 2.6 



2ML_,_,, , ,,/M|^ 



tan/? + cot a = — — f tan /? cos 2/? + C — 1|!- . (4.8) 



Mi ' \MX 



Thus, the first term in eq. ^^ and in eq. [4.7| is of order M'^y^ ian (3 / M\, 
and therefore decouples in the limit of large Ma- However, the approach to 
decoupling is delayed in the large tan/3 regime.^ Specifically, for values of 
M\ ~ Mz tan/?, we see that tan/3+cot a ~ 0{1). For example, if tan /3 ~ 50, 
then even for values of Ma ~ 1 TeV, decoupling has not yet set in. 



Other terms in eqs. ^]5| and ^^ also exhibit delayed decoupling. In par- 
ticular, eq. ^^ implies that 

n = X. + o(^^^25lgHil^). (4.9) 

SO that Yfe is also enhanced at large tan/3. Hence, all terms in eqs. ^^ and 



0| that are proportional to either X^ or Yf, are of order M'^y^ tan 13 / M^jj gy ■ 



Again, if tan/? ~ 50 and MgusY ~ 1 TeV, decoupling has not yet set in. 

The remaining terms in eqs. |4.5| and [4.7| exhibit the expected decoupling 
in the usual sense (with no delay). In particular, we may set a = /3 — 7r/2 in 
the decoupling limit to obtain 

cos/3sin(c. + /?) ^^^^ /ilg^\ 

sma y M4 / 

which exhibits no tan/? enhancement. All remaining factors of tan/? are 
multiplied by the appropriate power of m^, and since mfetan/3 ~ Mew^ no 
delayed decoupling results from these terms. 

We have thus shown analytically that the one-loop SQCD corrections to 
the h^hb coupling decouple in the limit of large Msusy and large Ma- The 
decoupling takes the generic form: 

AsQCD-C,^^ + C,j^^- (4.11) 

In general Ci approaches a non-zero constant as Msusy —* 00, while C2 ap- 
proaches a (different) non-zero constant as Ma -^ 00. Thus, the decoupling 
limit requires both Ma and Msusy to become simultaneously large (as com- 
pared to Mew)- However, we will demonstrate that in some cases the SQCD 
radiative corrections vanish in the limit where some SUSY particle masses 
are large, independent of the value of Ma- 

This decoupling is shown numerically in figs. ^ and |^. In fig. |^, we plot 

''The enhancement of the radiatively-corrected h'^bb coupHng at large tan/3 has also 
been emphasized in refs. P,[^,^p^|l5|. 

^In our numerical analysis we have taken the 6-quark pole-mass to be 4.75 GeV and 
as = 0.119. Because of the experimental constraints on the sbottom masses, we consider 
only regions of parameter space in which both sbottoms are heavier than 100 GeV. 
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Figure 2: Asqcd as a function of particle mass for tan/3 = 50 and Msusy = 
Ml = Mr = Ms = Mg = /i = Ab. The curves (a) are plotted vs. Msusy, 
with Ma = 3000 GeV; whereas the curves (b) are plotted vs. Ma, with 
Msusy = 3000 GeV. Solid lines are based on the exact one-loop formula and 
dashed lines are based on the analytic approximation of eq. [4.5|. 



the exact one-loop expression for Asqcd (solid lines) and the expansion of 
eq. |]5| (dashed lines) for tan/5 = 50 and Msusy = Ml = Mr = Ms = 
Mg = n = Aij. The hues labeled (a) show Asqcd as a function of Msusy- 
We have fixed Ma very large, Ma = 3000 GeV, in order to eliminate the 
contribution to Asqcd that decouples at large Ma- We use the exact tree- 
level formula for cot a as a function of Ma and tan f3. The lines labeled (b) 
show Asqcd as a function of Ma, where now we have fixed Msusy to be 
very large, Msusy = 3000 GeV, in order to examine only the contribution 
to Asqcd that does not decouple at large Msusy- We note that for very 
large Msusy and Ma = I TeV, Asqcd is of order -1% for tan/5 = 50. We 
have plotted our results for ^Mg positive. In the approximation of neglecting 
terms not enhanced by large tan/3 or cot a, changing the sign of ^Mg simply 
flips the sign of Asqcd- 

In fig. H we again plot the exact one-loop expression for Asqcd (solid 
lines) and the expansion of eq. ^]5| (dashed lines) for all the SUSY mass 
parameters equal, Msusy = Ml = Mr = Ms = Mg = /x = A^,, and three 
values of tan/?. Note the change in the vertical scale for the plots with 
different values of tan/?. We show the dependence of Asqcd on Msusy 
(left-hand panels) and Ma (right-hand panels). Clearly, in the limit of large 
Msusy, '^sqcd tends to a non- vanishing constant, and this constant tends 
to zero in the large Ma limit. Similarly, in the limit of large Ma, Asqcd 

^Although we have chosen Ml — Mn for simphcity, our results are not particularly 
sensitive to this choice as long as |M| — M^| ^ mfj\Xb\ {c.f. the remarks below eq. 4.5). 
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Figure 3: Asqcd for Msusy = Ml = Mr = Ms = Mg = fi = A^, with 
tan/3 = 8 (top panels), 30 (middle panels), and 50 (bottom panels). The 
solid lines are based on the exact one-loop formula and the dashed lines are 
based on the analytic approximation of eq. |4.5| . In the left-hand panels we 
plot Asqcd as a function of Msusy for Ma = 200, 300, and 500 GeV. In 
the right-hand panels we plot Asqcd as a function of Ma for Msusy = 200, 
300, and 500 GeV. For tan/? = 50, the value of Msusy = 200 GeV yields a 
negative mass-squared for the lighter sbottom, so this value is not shown in 
the bottom right panel. 



15 



tends to a non- vanishing constant, and this constant tends to zero in the 
large Msusy hmit. 

Notice that from the numerical comparison between the exact and ana- 
lytic formulae in fig. ||, we can conclude that our expansion is a good ap- 
proximation for large enough SUSY mass parameters. In particular, it is 
reasonably accurate for Msusy larger than 300 GeV. Also, it is clear that as 
tan P grows, not only does Asqcd increase in magnitude, but the agreement 
between the exact and analytic formulae becomes worse at low Msusy- This 
is due to the fact that the splitting between the squared masses of the two 
sbottoms in the maximal mixing case is proportional to mbtan/3, which we 
have taken to be of order Mew iii our expansion. As tan/5 increases, the 
mass of the lighter sbottom decreases, and the higher order terms that we 
have neglected in our expansion become more important. 

All numerical results presented so far correspond to fiMg > 0. In the case 
of fiMg < 0, the qualitative features of |A5qc'd| remain unchanged. From 
the analytic formulae derived in this section, one can see that at large tan/5 
the dominant effect of changing the sign of fiMg is to change the overall sign 
of Asqcd- We can illustrate this point in the simple limiting case in which 



all SUSY mass parameters and Ma are equal. Simplifying eq. [4.5| in this 
limit, we end up with a simple formula for the case of fiMg > 0: 

Asqcd = ^ {-|p— cos 2/5 (7tan/5 - 2) + ff (tan/5 - 1) 

6TT {OlVlsuSY ^-^^^^SUSY 

+||^('--3-4) + 0(^)|. (4.12) 

'^^^^SUSY \ ^^^SUSY / ) 

where Msusy = Ms = Mg = fi = Ai, = Ma- To obtain the result for 
fiMg < 0, one replaces tan/5 with — tan/5 in eq. |4.12|. The formula of eq. [4.12 



is plotted in fig. |^ for three values of tan /5 and both signs of ^ (taking Mg to 
be positive, by convention). Clearly, Asqcd decouples like {Mew/Msusy)j 
but this decoupling is delayed at large tan/5. For example, even at Msusy = 
1 TeV, IA5QC15I ^ 1% for tan/5 ~ 30. Note that as expected, changing the 
sign of /i simply changes the sign of the dominant contribution to Asqcd- 
In the remainder of our analysis, we will display results only for /x > 0. 

Next, we consider the decoupling of the SQCD corrections to the hPbb 
coupling as individual SUSY particles become heavy compared to the com- 
mon SUSY mass scale. We examine three cases: large Ms with maximal 
sbottom mixing, large Mg with maximal sbottom mixing, and one heavy 
sbottom state with near-zero sbottom mixing. 

We first consider the case of large Ms with maximal sbottom mixing, 
with Ms ^ Mg ~ yU ~ Aft :» Mew- If Ms is taken large while the rest of 
the SUSY mass parameters remain fixed, then we may expand the functions 
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Figure 4: Asqcd as a function of Msusy for Msusy = Ml = Mr = Ms = 
Mg = |/i| = Ah = Ma and tan/3 = 8, 30, 50. Both positive and negative fi 
cases are shown. The sohd hues are based on the exact one-loop formula and 
the dashed lines are based on the analytic approximation of eq. [4.12 . 



f^{R) in eq. U for M5 > Mg, or i? < 1. The result is: 
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, Mj COS p sm(a + p) ^. ,^ M^^ 

tan/3 + cot a) + — | ^-^ ^^ + ^ — I 

^ Ml sma ^ \M|, 



(4.13) 

where M is one of the lighter SUSY particle masses. Note that in this limit, 
the SQCD corrections decouple like M'^/Mg even for light M^. Thus it is 
only in the case of large Mg and /i, of the same order as Ms, that large M^ 
is required for decoupling. In fig. |^ we plot the exact one-loop expression 
for Asqcd and the expansions of eqs. ^^ and |4.13| as a function of Ms, for 



fixed Mg = ^ = Ab = Ma = 200 GeV and three different values of tan/?. 
This figure shows the decoupling of Asqcd with Ms as discussed above. 

Similarly we examine the case of a very heavy gluino compared to the 
rest of the SUSY spectrum. We still focus on the case of maximal sbottom 
mixing. Expanding the functions fi{R) in eq. |4.5| for Mg ^ Ms, or i? ^ 1, 
we see that in this case the SQCD corrections decouple with the gluino mass 
like M/Mg, where again M is one of the other light SUSY masses: 
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(4.14) 



M2 



Note that the decoupling of the SQCD corrections at large Mg (with all other 
SUSY mass parameters held fixed) is very slow: Asqcd falls off only as one 



power of Mg. This is due to the factor of Mg in the numerator of eqs. |3.6| and 
p.8| , which arises from the gluino propagator. Asqcd is also enhanced by the 
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Figure 5: Asqcd as a function of Ms (assuming Ml = Mr = Ms) for 
Mg = fi = Af, = Ma = 200 GeV and tan/3 = 8, 30, 50 (top to bottom). Solid 
lines are based on the exact one-loop expression, dashed lines are based on 
the analytic expansion of eq. Ol, and dotted lines are based on the large Ms 



expansion of eq. [4.13 . 
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Figure 6: Asqcd as a function of Mg for M^ = Mr = Ms = ^ = Ah = Ma = 
300 GeV and tan/3 = 8,30,50. Solid lines are based on the exact one-loop 
expression, dashed lines based on are the analytic expansion of eq. |4.5|, and 



dotted lines are based on the expansion for large Ms of eq. 4.14. 
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factor log(M|/M|). We again see the phenomenon of delayed decouphng at 
large tan/? due to the terms in eq. f4.14| proportional to either Xb or Yt. 

In fig. ^ we plot the exact one-loop expression for Asqcd and the ex- 
pansions of eqs. [4.5| and |4.14| as a function of Mg, for Ms = n = Ah = 
Ma = 300 GeV and three different values of tan/3. This figure shows the 
slow decoupling of Asqcd with Mg. For example, for Mg = 500 GeV and 
tan/5 = 30, Asqcd ^ -6% for Ms = fi = Ab = Ma = 300 GeV. If the 
latter masses are reduced to 200 GeV, one finds Asqcd — —13%, which is a 
significant correction. Fig. ^ also illustrates the validity of the large gluino 
mass expansion. This expansion is particularly poor for large values of tan /? 
out to a very large gluino mass of about 2000 GeV. 

Finally we study the case in which one of the sbottoms becomes heavy 
while the other sbottom mass and the rest of the SUSY mass parameters are 
fixed. We choose Mr-:^ Ml -^ Mg ^ ^i 



Ab > Mew, so that M^^ > M^^. 



This is necessarily the case of near-zero sbottom mixing. Expanding eq. ^J_ 
in inverse powers of M^ , we find: 
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b2 



b2. 



where again M is one of the other light SUSY masses and the function h{Ri) 
is given in the Appendix. Note that the first term does not decouple as 
Ml is taken large. This behavior is independent of the value of Ma (and 
therefore holds even if Ma — > oo). However, this first term is not enhanced 
by large tan/3 and is numerically negligible as can be seen in fig. |^. The 
terms that are enhanced by large tan/? decouple like M'^/M^ . In fig. |^ we 



plot the exact one-loop expression for Asqcd and the expansions of eqs. ^ 
and ^J^, as a function of M^^, for Mj,^ = Mg = fx = Ab = Ma = 200 
GeV and three different values of tan/3. Clearly, in order for Asqcd to be 
large in the case of near- zero sbottom mixing, both of the sbottoms must 
be reasonably light. Note however that, due to the enhancement in tan/?, 
the 1/M^ suppression is not so small. As an example, for tan/3 = 50, 
Mi^ =mI = ii = Ab = MA = 200 GeV and M^^ = 500 GeV [1000 GeV], one 
obtains Asqcd ^ -10% [-5%]. 

The various cases examined in this section can be summarized by spec- 
ifying the behavior of Ci and C2 of eq. |4.11| on the model parameters. In 

is identified with the 



Table 1, four cases are shown. In all cases, M- 



SUSY 
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Figure 7: Asqcd as a function of M^^ , with M^^ = Mg = ft = Af, = Ma = 
200 GeV and tan/5 = 8,30,50. Solid lines are based on the exact one- 
loop formula, dashed lines based on are the analytic expansion for near-zero 



sbottom mixing of eq. 4.7, and dotted lines are based on the expansion for 



large M^^ of eq. |4.15 



largest supersymmetry-breaking mass, while M refers to a possible interme- 
diate supersymmetric mass scale (with Mew <^ M <^ Msusy)- The presence 
of a factor of tan/5 (unless multiplied by M/ Msusy) indicates delayed de- 
coupling. In the case of Mg = Msusy, C2 ~ {Msusy /M) tan P implies a 

Fi- 



delayed decoupling that vanishes only as a single power of 1/Msusy- 
nally, in the case of large M^^, C2 ~ Msusy/ ^^ implies no decoupling as 
Msusy -^ 00 with M held fixed. This is not a violation of the usual decou- 
pling theorem 0,0], since in the latter case, only part of the supersymmetric 
spectrum has been removed from the low-energy effective theory. Decoupling 
is recovered in the limit of M — ;► 00, as expected. 



Case 



b mixing 



Ci 



C, 



Ms^Mg = Msusy 
Ms = Msusy > M 
Mg = Msusy > M 
Mi^ = Msusy > M 



maximal 
maximal 
maximal 
near-zero 



tan/5 



(MVM|^5^) tan/? 



tan/5 
1 

{M/ Msusy) tan /5 (Msusy /M) tan /5 
(MVM|^5^) tan /5 MI^sy/M' 



Table 1: Approach to decoupling of the one-loop 0{as) radia- 
tive corrections to the h'^bb vertex: Asqcd ~ C*i(M|;,y/M^) + 
C2{M^-^^/ Msusy)- See text for further discussion. 
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5 Conclusions 

In this paper we have studied the one loop SQCD corrections to the h^bb cou- 
phng in the hmit of large SUSY masses. In order to understand analytically 
the behavior of the corrections in this limit, we have performed expansions 
for the SUSY mass parameters large compared to the electroweak scale. We 
have shown that for the SUSY mass parameters and Ma large and all of the 
same order, the SQCD corrections decouple like the inverse square of these 
mass parameters. However, if the mass parameters are not all of the same 
size, then this behavior can be modified. If Ma is light, then the SQCD 
corrections to the h^bb coupling generically do not decouple in the limit of 
large SUSY mass parameters. In this case, the low-energy theory at the elec- 
troweak scale contains two full Higgs doublets with Higgs-fermion couplings 
of the general type-Ill model. 

We have also examined three cases in which there is a hierarchy among 
the SUSY mass parameters. In the case of maximal sbottom mixing with 
Ms large and the other SUSY mass parameters and Ma of order a common 
mass scale M (chosen such that Mew <^ M -C Ms), the SQCD corrections 
decouple like M'^/Mg. Second, we examined the case of a large gluino mass 
with the other SUSY mass parameters of order a common mass scale M 
(chosen such that Mew <^ M <^ ^g)- ^^ this case we found that the 
SQCD corrections decouple more slowly, like (M/M^) log (M|/M|). Finally, 
we examined the case in which one sbottom is much heavier than the other 
SUSY mass parameters, which are fixed at a scale M. In this case the mixing 
angle in the sbottom sector is near zero. We found that the piece of the SQCD 
corrections that is enhanced at large tan/3 decouples like M'^/M^ . There is 
also a piece of the SQCD corrections that does not decouple as M^^ is taken 
large, but it is not enhanced by tan/3 and is numerically negligible compared 
to the decoupling piece, up to a very high value of the heavier sbottom mass. 

The decoupling behavior of the SQCD corrections to the h^bb coupling im- 
plies that distinguishing the lightest MSSM Higgs boson from the SM Higgs 
boson will be very difficult if A^ and the SUSY spectrum are heavy, even 
after one-loop SUSY corrections are taken into account. However, because 
of the enhancement at large tan/3, the onset of decoupling is delayed, and 
the corrections can still be at the percent level for tan/3 ~ 50 and all SUSY 
mass parameters and Ma of order 1 TeV. If one or both of the sbottoms, the 
gluino, and/or A^ lie below the TeV scale, then the SQCD corrections will 
be larger still. The decoupling limit provides a challenge for Higgs searches 
at future colliders. Even if the light CP-even Higgs boson is found, the di- 
rect discovery of supersymmetric particles may be essential for unraveling 
the origin of electroweak symmetry breaking. 
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Appendix 
A Expansions of loop functions 

In this Appendix we define our notation for the two- and three-point integrals 



that appear in eqs. |3.6| and |3.8| and give formulae for their expansions in 



powers of the SUSY mass parameters. 

We follow the definitions and conventions of |34]. The two-point integrals 
are given by: 

>^'-" I (gsF^T^strip^ - li^ («»^ '"^'' <'^^ "'• '"^' ■ 

(A.l) 
The derivatives of the two-point functions are defined as follows: 

B'.^^ip'-mlml) = ^Bo,i{q';mlml) . (A.2) 

Finally, the three-point integrals are given by: 



/x^-^ 



d^k {l;k^} 



(27r)^ [P - mj] [{k + piY - ml] [{k + pi + pa)' - ml] 

{Co;p^Cii+pi^Ci2}{pl,pl,p'^;ml,ml,ml), (A.3) 



167r2 

where p = — pi — p2 ■ 

We now give the large Msusy expansions of the loop integrals. 
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A.l Maximal bi — bji mixing 

The loop integrals are expanded as follows, where M^ and A^ are defined in 



eq. |4.4 



Co{mlMlo,mlMlMl,Ml: 



1 A^ A^ M^ 



21^^-"' ' SMl'^'^"' AMl'-'"-^^' 24M| 



Co{mlMl,,ml-MlMl,Ml 



9' "'"b2' te'^ 
1 „ , „. A2 



r/i(i?) 



Co{mlMlo,ml,MlMl,Ml) 
1 A"^ 

" ~2M|^'^^^ ~ T2M| 



/.(i?) - ^/3(i?) 









CuimlMlo,mlMlMl,Ml) 



/7(i?) 



Cu{mlMlo,mlMlMl,Ml 
1 A2 



3M| 



A4 /M-2 



B,{ml-MlMl)-B,{ml-MlMl)c^ 



B'.iml Ml Ml) - B',{ml Ml Ml: 



A^ 



A2 



5U<; M^, M^J + 5Umf ; M^, M^ 



b2^ 



6M| 



6M| 



/8(i?) 



/4(i?) 



A^ 



15M| 



/9(i?). (A.4) 



The functions fi{R) are given in terms of the ratio R = Mg/Ms'. 
2 



MR) 
MR) 
MR) 
MR) 



;i-i?2)2 

3 

;i-i?2)3 

4 

;i-i?2)4 

4 

;i - i?2)4 

3 

;i-i?2)3 

5 

;i-/?2)5 

5 



l-i?' + i?'logi?' 



l-i?^ + 2i?'logi?' 



1 + f /?" - 3i?* + i/?'' + 3R' log i?' 



i-3i?2 + |/?4^i?6_3^4iQg^2 



l-2R^ + lR^-RUogR' 



I - 4i?2 + 4i?6 - |i?^ - 6i?^ log i?2 



i - 2i?^ + 6i?" - f i?^ - i?« + 4i?'' log R' 
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f9{R) 



12 



(l-i?2)4 

5 

(1 - /?2)6 



i + 2i?2 _ 5^4 ^ 2i?2 log ^2 ^ ^4 iQg ^2 



1 - UR-" - 36i?^ + 44i?'' + 3i?^ 



-24i?^logi?^-36i?^logi?^ 



(A.5) 



Note that in the special case Mg = Ms, R = 1 and /«(!) = 1. 



A. 2 Near-zero bi — bji mixing 

The loop integrals are expanded as follows: 

1 



CoK,M,^o,m^;M^,M(,M^; 



Co{mlM',o,ml,MlMl,Ml] 



CoimlMlo,ml;MlMl,Ml] 



Cu{mlMlo,ml;MlMl,Ml 



CuimlMlo,ml;MlMl,Ml: 



2M-^ 
1 



2 /l(^l) - o..'!r4 /4(^1^ 



2M^ 



/l(^2 



24M^ 
24Mi^ 

^2 



/4(i?2) 



hiR,,R2) ^Mloh2iRi,R2) 






3M-, 



3Mi 
62 



30Mi 

bi 

/5(^2) + i^KVTrfriRi 



30M^ 

^2 



BoK^; M|, MfJ -5oK;M|,Mj) ^ -/ii(i?i,i?2) 



B',{mt;M^^,Mi^)-B',{mi;M^,Mi 



1 



fe^ 6M^ 



■/2(i?i) 



6M? 

b2 



/2(i?2) 



B[imi;M^^,Ml) + B[imt;M^,Ml 



l2M-f 



-h{Ri) 



12M1 
02 



/4(i?2),(A.6) 



where Ri = Mg/M^, {i = 1,2). The functions fi{R) were given in eq. |A.5 
The functions hi{Ri,R2) and h2{Ri,R2) are defined as follows: 



hi{Ri, R2) = h{Ri) - /i(i?2), with h{R) 



\ogR^ 
l-i?2' 



Mi?i,i?2) = l + 5^^^4^;^ 



2d2 

2 



2(l-i??)(l-i?i 



logi?? 



2(i?? - Rl) 



[l-R^.y 



[Rl + -R2 — '^Ri) 



_\ogR? 



{Rl + R2- 2R2) 



(A.7) 
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The functions h and /i2 have the following properties: 
h2{Ri,R2) = h2{R2,Ri), 



hoihRo 



0--R 



2\2 



I- Rl- lRl+ (\ + RD^ogRl] . (A.8) 



References 

[1] H. E. Haber and Y. Nir, Nucl. Phys. B335, 363 (1990); H. E. Haber, 
in Proceedings of the US-Polish Workshop, Warsaw, Poland, Septem- 
ber 21-24, 1994, edited by P. Nath, T. Taylor, and S. Pokorski (World 
Scientific, Singapore, 1995) pp. 49-63. 

[2] A. Dobado, M. J. Herrero and S. Penaranda, Eur. Phys. J. C7, 313 
(1999); Eur. Phys. J. C12, 673 (2000); |E^2^£E70002T34[ 



[3] 
[4] 

[5] 
[6] 

[7 

[8 

[9 

[10 

[11 

[12 
[13 



T. Appelquist and J. Carazzone, Phys. Rev. Dll, 2856 (1975). 

M. Carena, S. Mrenna and C. E. M. Wagner, Phys. Rev. D60, 075010 
(1999). 



M. Carena, S. Mrenna and C. E. M. Wagner, |hep-ph/9907422 



D. Rainwater, D. Zeppenfeld and K. Hagiwara, Phys. Rev. D59, 014037 
(1999); T. Plehn, D. Rainwater and D. Zeppenfeld, Phys. Rev. D61, 
093005 (2000). 

F. Borzumati, G. R. Farrar, N. Polonsky and S. Thomas, Nucl. Phys. 
B555, 53 (1999). 

K.S. Babu and C. Kolda, Phys. Lett. B451, 77 (1999). 

A. Dabelstein, Nucl. Phys. B456, 25 (1995). 

J. A. Coarasa, R. A. Jimenez and J. Sola, Phys. Lett. B389, 312 (1996). 

L. J. Hall, R. Rattazzi and U. Sarid, Phys. Rev. D50, 7048 (1994); 
R. Hempfling, Phys. Rev. D49, 6168 (1994); M. Carena, M. Olechowski, 
S. Pokorski and C. E. M. Wagner, Nucl. Phys. B426, 269 (1994). 

D. M. Pierce, J. A. Bagger, K. Matchev and R.-J. Zhang, Nucl. Phys. 
B491, 3 (1997). 



S. Heinemeyer, W. Hollik and G. Weiglein, |hep-ph/0003"022 . 



25 



[14] G.L. Kane, G.D. Kribs, S.P. Martin and J.D. Wells, Pliys. Rev. D53, 
213 (1996); H. Baer and J.D. Wells, Pliys. Rev. D57, 4446 (1998). 

[15] W. Loinaz and J.D. Wells, Phys. Lett. B445, 178 (1998). 

[16] E. Braaten and J. P. Leveille, Phys. Rev. D22, 715 (1980); N. Sakai, 
Phys. Rev. D22, 2220 (1980); T. Inami and T. Kubota, Nucl. Phys. 
B179, 171 (1981). 

[17] M. Carena, D. Garcia, U. Nierste and C. E. M. Wagner, Nucl. Phys. 
B577, 88 (2000). 



[18] H. Eberl, K. Hidaka, S. Kraml, W. Majerotto and Y. Yamada, |hep- 
I ph/9912463| 



[19] J. F. Gunion, H. E. Haber, G. Kane and S. Dawson, The Higgs Hunter's 
Guide (Addison-Wesley, Reading, MA, 1990) [errata: |hep-ph/ 9302772 



[20] J. F. Gunion and H. E. Haber, Nucl. Phys. B272, 1 (1986); B278, 449 
(1986) [E: B402, 567 (1993)]. 

[21] H. E. Haber and R. Hempfling, Phys. Rev. Lett. 66, 1815 (1991); Phys. 
Rev. D48, 4280 (1993). Y. Okada, M. Yamaguchi and T. Yanagida, 
Prog. Theor. Phys. 85, 1 (1991); Phys. Lett. B262, 54 (1991); J. Ellis, 
G. Ridolfi and F. Zwirner, Phys. Lett. B257, 83 (1991); Phys. Lett. 
B262, 477 (1991); R. Barbieri and M. Frigeni, Phys. Lett. B258, 167 
(1991); Phys. Lett. B258, 395 (1991); 

[22] R. Hempfling and A. H. Hoang, Phys. Lett. B331, 99 (1994); 
J. A. Casas, J. R. Espinosa, M. Quiros and A. Riotto, Nucl. Phys. B436, 
3 (1995) [E: B439, 466 (1995)]; M. Carena, J. R. Espinosa, M. Quiros 
and C. E. M. Wagner, Phys. Lett. B355, 209 (1995); M. Carena, 
M. Quiros and C. E. M. Wagner, Nucl. Phys. B461, 407 (1996); 
H. E. Haber, R. Hempfling and A. H. Hoang, Z. Phys. C75, 539 (1997). 

[23] J. R. Espinosa and R.-J. Zhang, JHEP 0003, 026 (2000); ^eg 
I ph/00032lB| . 

[24] P. H. Chankowski, S. Pokorski and J. Rosiek, Phys. Lett. B274, 191 
(1992); Nucl. Phys. B423, 437 (1994); A. Dabelstein, Z. Phys. C67, 

495 (1995). 

[25] S. Heinemeyer, W. Hollik and G. Weiglein, Phys. Rev. D58, 091701 
(1998); Phys. Lett. B440, 296 (1998); Eur. Phys. J. C9, 343 (1999); 
Phys. Lett. B455, 179 (1999). 

M. Carena, H. E. Haber, S. Heinemeyer, W. Hollik, C. E. M. Wagner 
and G. Weiglein, Nucl. Phys. B580, 29 (2000). 



26 



[27] P. Bock et al. [ALEPH, DELPHI, L3 and OPAL Collaborations and 
the LEP working group for Higgs boson searches], CERN-EP-2000-055 
(April, 2000). 

[28] T. Affolder et al. [CDF Collaboration], Phys. Rev. Lett. 84, 5704 (2000). 

[29] The ALEPH Collaboration, ALEPH 2000-012 CONF 2000-009, submit- 
ted to the 2000 Winter Conferences. 

[30] S. Abachi et al. [DO Collaboration], Phys. Rev. Lett. 75, 618 (1995); F. 
Abe et al. [CDF Collaboration], Phys. Rev. D56, R1357 (1997). 

[31] D. Atwood, L. Reina and A. Soni, Phys. Rev. D55, 3156 (1997) and 
references therein. 

[32] L. Randall and R. Sundrum, Nucl. Phys. B557, 79 (1999); G. F. Giu- 
dice, M. A. Luty, H. Murayama, and R. Rattazzi, JHEP 9812, 027 
(1998). 

[33] G. D. Kribs, Phys. Rev. D62, 015008 (2000). 

[34] W. Hollik, in Precision Tests of the Standard Electroweak Model, edited 
by P. Langacker (World Scientific, Singapore, 1995), p. 37-116. 



27 



